This is the analysis script to prepare the data to run analyses for the project “Anxiety, depression and trauma symptom change during the COVID-19 pandemic: retrospective versus objective assessment” - Young et al (2020)
Script written by K Purves, K Thompson, C Huebel and M Davies. Email: kirstin.purves@kcl.ac.uk, katie.thompson@kcl.ac.uk, christopher.1.huebel@kcl.ac.uk, molly.davies@kcl.ac.uk
#Set up
Clear global environment
Retrieve the current date to use it for file endings to not overwrite files when one exports files
#Packages Install packages (if they are not available in your version of R)
Load packages
Registered S3 method overwritten by 'pryr':
method from
print.bytes Rcpp
For best results, restart R session and update pander using devtools:: or remotes::install_github('rapporter/pander')
Attaching package: 'polycor'
The following object is masked from 'package:psych':
polyserial
corrplot 0.84 loaded
── Attaching packages ─────────────────────────────────────── tidyverse 1.3.0 ──
✓ ggplot2 3.3.2 ✓ purrr 0.3.4
✓ tibble 3.0.4 ✓ dplyr 1.0.2
✓ tidyr 1.1.2 ✓ stringr 1.4.0
✓ readr 1.3.1 ✓ forcats 0.5.0
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
x ggplot2::%+%() masks psych::%+%()
x ggplot2::alpha() masks psych::alpha()
x dplyr::filter() masks stats::filter()
x dplyr::lag() masks stats::lag()
x tibble::view() masks summarytools::view()
Define colours for plotting this are the standard coping colours
Choose in this chunk which palette to use 04.07.2020 - Default to 2 colour COPING palette
Set up ggplot theme for the plots
[1] 65223 643
[1] "gad.sum_score_prepan" "pcl.sum_score_prepan"
[3] "phq.sum_score_prepan" "phq.sum_score_8items_prepan"
[5] "ocir.sum_score_prepan" "gad.sum_score_base"
[7] "gad.sum_score_retro" "pcl.sum_score_base"
[9] "pcl.sum_score_retro" "phq.sum_score_base"
[11] "phq.sum_score_8items_base" "phq.sum_score_retro"
[13] "phq.sum_score_8items_retro" "ocir.sum_score_base"
[15] "ocir.sum_score_retro"
[1] "phq.sum_score_prepan" "phq.sum_score_base" "phq.sum_score_retro"
[1] "gad.sum_score_prepan" "gad.sum_score_base" "gad.sum_score_retro"
[1] "ocir.sum_score_prepan" "ocir.sum_score_base" "ocir.sum_score_retro"
[1] "pcl.sum_score_prepan" "pcl.sum_score_base" "pcl.sum_score_retro"
[1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan"
[3] "phq.diff_score_base_retro" "gad.diff_score_retro_prepan"
[5] "gad.diff_score_base_prepan" "gad.diff_score_base_retro"
[7] "pcl.diff_score_base_prepan" "ocir.diff_score_retro_prepan"
[9] "ocir.diff_score_base_prepan" "ocir.diff_score_base_retro"
[1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan"
[3] "phq.diff_score_base_retro"
[1] "gad.diff_score_retro_prepan" "gad.diff_score_base_prepan"
[3] "gad.diff_score_base_retro"
[1] "ocir.diff_score_retro_prepan" "ocir.diff_score_base_prepan"
[3] "ocir.diff_score_base_retro"
Exclude GLAD/EDGI individuals that didn’t take part in COPING Exclude individuals with NA for the MHD categories Individuals with NAs for age/gender have already been excluded in pre-processing script, but they get excluded here to be safe and to double check previous programming.
[1] 65223 643
[1] 34465 643
[1] 30758
Frequencies
dat$gad.pandemic_felt_feelings
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------ --------- -------------- --------- --------------
Much worse 10 15.38 15.38 15.38 15.38
A little worse 32 49.23 64.62 49.23 64.62
No different 16 24.62 89.23 24.62 89.23
A little better 5 7.69 96.92 7.69 96.92
Much better 2 3.08 100.00 3.08 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 2474 19.56 19.56 19.55 19.55
A little worse 5321 42.08 61.64 42.05 61.61
No different 3599 28.46 90.10 28.44 90.05
A little better 854 6.75 96.85 6.75 96.80
Much better 398 3.15 100.00 3.15 99.94
<NA> 7 0.06 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 879 6.72 6.72 6.71 6.71
A little worse 3853 29.45 36.17 29.42 36.13
No different 7897 60.37 96.54 60.30 96.43
A little better 336 2.57 99.11 2.57 99.00
Much better 117 0.89 100.00 0.89 99.89
<NA> 14 0.11 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------ --------- -------------- --------- --------------
Much worse 2357 27.35 27.35 27.25 27.25
A little worse 3568 41.41 68.76 41.24 68.49
No different 2128 24.70 93.45 24.60 93.09
A little better 363 4.21 97.67 4.20 97.28
Much better 201 2.33 100.00 2.32 99.61
<NA> 34 0.39 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$phq.pandemic_felt_feelings
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------ --------- -------------- --------- --------------
Much worse 12 18.46 18.46 18.46 18.46
A little worse 27 41.54 60.00 41.54 60.00
No different 17 26.15 86.15 26.15 86.15
A little better 6 9.23 95.38 9.23 95.38
Much better 3 4.62 100.00 4.62 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 2282 18.09 18.09 18.04 18.04
A little worse 5645 44.74 62.83 44.61 62.65
No different 3217 25.50 88.33 25.42 88.07
A little better 1034 8.20 96.53 8.17 96.25
Much better 438 3.47 100.00 3.46 99.71
<NA> 37 0.29 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
Much worse 899 6.90 6.90 6.86 6.86
A little worse 4305 33.04 39.94 32.87 39.74
No different 7351 56.42 96.37 56.13 95.87
A little better 338 2.59 98.96 2.58 98.45
Much better 135 1.04 100.00 1.03 99.48
<NA> 68 0.52 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------- ------ --------- -------------- --------- --------------
Much worse 2288 26.61 26.61 26.45 26.45
A little worse 3768 43.83 70.44 43.56 70.00
No different 1885 21.93 92.37 21.79 91.79
A little better 421 4.90 97.27 4.87 96.66
Much better 235 2.73 100.00 2.72 99.38
<NA> 54 0.62 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$Ethnicity_collapsed
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------ --------- -------------- --------- --------------
European 60 92.31 92.31 92.31 92.31
Mixed or multiple ethnic origins 3 4.62 96.92 4.62 96.92
Asian or Asian British 1 1.54 98.46 1.54 98.46
African or African British 0 0.00 98.46 0.00 98.46
Other 1 1.54 100.00 1.54 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
European 12090 95.78 95.78 95.55 95.55
Mixed or multiple ethnic origins 276 2.19 97.96 2.18 97.73
Asian or Asian British 98 0.78 98.74 0.77 98.51
African or African British 37 0.29 99.03 0.29 98.80
Other 122 0.97 100.00 0.96 99.76
<NA> 30 0.24 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
European 9927 97.16 97.16 75.80 75.80
Mixed or multiple ethnic origins 118 1.15 98.32 0.90 76.70
Asian or Asian British 126 1.23 99.55 0.96 77.66
African or African British 46 0.45 100.00 0.35 78.02
Other 0 0.00 100.00 0.00 78.02
<NA> 2879 21.98 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------ --------- -------------- --------- --------------
European 8124 94.32 94.32 93.91 93.91
Mixed or multiple ethnic origins 174 2.02 96.34 2.01 95.92
Asian or Asian British 169 1.96 98.30 1.95 97.87
African or African British 76 0.88 99.19 0.88 98.75
Other 70 0.81 100.00 0.81 99.56
<NA> 38 0.44 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$Gender_collapsed
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------- ------ --------- -------------- --------- --------------
Male 1 1.54 1.54 1.54 1.54
Female 62 95.38 96.92 95.38 96.92
Non-binary/Self-defined 2 3.08 100.00 3.08 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------- ------- --------- -------------- --------- --------------
Male 2279 18.01 18.01 18.01 18.01
Female 10127 80.04 98.05 80.04 98.05
Non-binary/Self-defined 247 1.95 100.00 1.95 100.00
<NA> 0 0.00 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------- ------- --------- -------------- --------- --------------
Male 5641 43.07 43.07 43.07 43.07
Female 7416 56.63 99.70 56.63 99.70
Non-binary/Self-defined 39 0.30 100.00 0.30 100.00
<NA> 0 0.00 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------- ------ --------- -------------- --------- --------------
Male 1777 20.54 20.54 20.54 20.54
Female 6782 78.40 98.94 78.40 98.94
Non-binary/Self-defined 92 1.06 100.00 1.06 100.00
<NA> 0 0.00 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$age_category_collapsed
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
16-18 7 10.77 10.77 10.77 10.77
19-25 17 26.15 36.92 26.15 36.92
26-35 21 32.31 69.23 32.31 69.23
36-45 8 12.31 81.54 12.31 81.54
46-55 4 6.15 87.69 6.15 87.69
56-65 5 7.69 95.38 7.69 95.38
66-70 3 4.62 100.00 4.62 100.00
71-75 0 0.00 100.00 0.00 100.00
76+ 0 0.00 100.00 0.00 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
16-18 341 2.70 2.70 2.70 2.70
19-25 1537 12.15 14.84 12.15 14.84
26-35 2850 22.52 37.37 22.52 37.37
36-45 2459 19.43 56.80 19.43 56.80
46-55 2821 22.30 79.10 22.30 79.10
56-65 1908 15.08 94.18 15.08 94.18
66-70 433 3.42 97.60 3.42 97.60
71-75 231 1.83 99.42 1.83 99.42
76+ 73 0.58 100.00 0.58 100.00
<NA> 0 0.00 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
16-18 10 0.08 0.08 0.08 0.08
19-25 291 2.22 2.30 2.22 2.30
26-35 1165 8.90 11.19 8.90 11.19
36-45 1440 11.00 22.19 11.00 22.19
46-55 2644 20.19 42.38 20.19 42.38
56-65 3901 29.79 72.17 29.79 72.17
66-70 1967 15.02 87.19 15.02 87.19
71-75 1265 9.66 96.85 9.66 96.85
76+ 413 3.15 100.00 3.15 100.00
<NA> 0 0.00 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
16-18 697 8.06 8.06 8.06 8.06
19-25 772 8.92 16.98 8.92 16.98
26-35 1076 12.44 29.42 12.44 29.42
36-45 874 10.10 39.52 10.10 39.52
46-55 1519 17.56 57.08 17.56 57.08
56-65 2201 25.44 82.52 25.44 82.52
66-70 786 9.09 91.61 9.09 91.61
71-75 514 5.94 97.55 5.94 97.55
76+ 212 2.45 100.00 2.45 100.00
<NA> 0 0.00 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$employment_prior_covid
Type: Character
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------ --------- -------------- --------- --------------
Employed 16 25.40 25.40 24.62 24.62
Key worker 23 36.51 61.90 35.38 60.00
Retired 5 7.94 69.84 7.69 67.69
Student 8 12.70 82.54 12.31 80.00
Unemployed 11 17.46 100.00 16.92 96.92
<NA> 2 3.08 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
Employed 3164 26.97 26.97 25.01 25.01
Key worker 5357 45.66 72.62 42.34 67.34
Retired 1435 12.23 84.85 11.34 78.68
Student 673 5.74 90.59 5.32 84.00
Unemployed 1104 9.41 100.00 8.73 92.73
<NA> 920 7.27 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
Employed 3398 26.47 26.47 25.95 25.95
Key worker 4414 34.38 60.85 33.70 59.65
Retired 4681 36.46 97.31 35.74 95.40
Student 106 0.83 98.14 0.81 96.20
Unemployed 239 1.86 100.00 1.82 98.03
<NA> 258 1.97 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------ --------- -------------- --------- --------------
Employed 2345 28.45 28.45 27.11 27.11
Key worker 2765 33.54 61.99 31.96 59.07
Retired 1913 23.21 85.20 22.11 81.18
Student 720 8.73 93.93 8.32 89.50
Unemployed 500 6.07 100.00 5.78 95.28
<NA> 408 4.72 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$anxiety_disorders
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------- ------ --------- -------------- --------- --------------
No anxiety disorder 27 41.54 41.54 41.54 41.54
Anxiety disorder 38 58.46 100.00 58.46 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------- ------- --------- -------------- --------- --------------
No anxiety disorder 1637 13.95 13.95 12.94 12.94
Anxiety disorder 10098 86.05 100.00 79.81 92.74
<NA> 918 7.26 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------- ------- --------- -------------- --------- --------------
No anxiety disorder 11303 86.49 86.49 86.31 86.31
Anxiety disorder 1765 13.51 100.00 13.48 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------- ------ --------- -------------- --------- --------------
No anxiety disorder 4210 55.69 55.69 48.66 48.66
Anxiety disorder 3350 44.31 100.00 38.72 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$depressive_disorders
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------------------- ------ --------- -------------- --------- --------------
No depressive disorder 18 27.69 27.69 27.69 27.69
Depressive disorder 47 72.31 100.00 72.31 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------------------- ------- --------- -------------- --------- --------------
No depressive disorder 1197 9.46 9.46 9.46 9.46
Depressive disorder 11455 90.54 100.00 90.53 99.99
<NA> 1 0.01 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------------------- ------- --------- -------------- --------- --------------
No depressive disorder 10312 78.91 78.91 78.74 78.74
Depressive disorder 2756 21.09 100.00 21.04 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------------------- ------ --------- -------------- --------- --------------
No depressive disorder 3734 49.39 49.39 43.16 43.16
Depressive disorder 3826 50.61 100.00 44.23 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$depression_and_anxiety
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------------- ------ --------- -------------- --------- --------------
No depressive or anxiety disorder 16 24.62 24.62 24.62 24.62
Depressive and anxiety disorder 36 55.38 80.00 55.38 80.00
Only depressive disorder 11 16.92 96.92 16.92 96.92
Only anxiety disorder 2 3.08 100.00 3.08 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------------- ------- --------- -------------- --------- --------------
No depressive or anxiety disorder 275 2.34 2.34 2.17 2.17
Depressive and anxiety disorder 9257 78.88 81.23 73.16 75.33
Only depressive disorder 1362 11.61 92.83 10.76 86.10
Only anxiety disorder 841 7.17 100.00 6.65 92.74
<NA> 918 7.26 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------------- ------- --------- -------------- --------- --------------
No depressive or anxiety disorder 9666 73.97 73.97 73.81 73.81
Depressive and anxiety disorder 1119 8.56 82.53 8.54 82.35
Only depressive disorder 1637 12.53 95.06 12.50 94.85
Only anxiety disorder 646 4.94 100.00 4.93 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------------- ------ --------- -------------- --------- --------------
No depressive or anxiety disorder 2935 38.82 38.82 33.93 33.93
Depressive and anxiety disorder 2551 33.74 72.57 29.49 63.41
Only depressive disorder 1275 16.87 89.43 14.74 78.15
Only anxiety disorder 799 10.57 100.00 9.24 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$eating_disorders
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------ ------ --------- -------------- --------- --------------
No eating disorder 4 6.15 6.15 6.15 6.15
Eating disorder 61 93.85 100.00 93.85 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------ ------- --------- -------------- --------- --------------
No eating disorder 11255 89.05 89.05 88.95 88.95
Eating disorder 1384 10.95 100.00 10.94 99.89
<NA> 14 0.11 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------ ------- --------- -------------- --------- --------------
No eating disorder 12897 98.51 98.51 98.48 98.48
Eating disorder 195 1.49 100.00 1.49 99.97
<NA> 4 0.03 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------ ------ --------- -------------- --------- --------------
No eating disorder 6592 93.19 93.19 76.20 76.20
Eating disorder 482 6.81 100.00 5.57 81.77
<NA> 1577 18.23 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$obsessive_compulsive_disorders
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------ --------- -------------- --------- --------------
No obsessive compulsive disorder 47 72.31 72.31 72.31 72.31
Obsessive compulsive disorder 18 27.69 100.00 27.69 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
No obsessive compulsive disorder 11025 87.14 87.14 87.13 87.13
Obsessive compulsive disorder 1627 12.86 100.00 12.86 99.99
<NA> 1 0.01 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
No obsessive compulsive disorder 12923 98.89 98.89 98.68 98.68
Obsessive compulsive disorder 145 1.11 100.00 1.11 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------- ------ --------- -------------- --------- --------------
No obsessive compulsive disorder 6953 91.97 91.97 80.37 80.37
Obsessive compulsive disorder 607 8.03 100.00 7.02 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$psychotic_disorders
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------- ------ --------- -------------- --------- --------------
No psychotic disorder 65 100.00 100.00 100.00 100.00
Psychotic disorder 0 0.00 100.00 0.00 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------- ------- --------- -------------- --------- --------------
No psychotic disorder 12238 96.83 96.83 96.72 96.72
Psychotic disorder 401 3.17 100.00 3.17 99.89
<NA> 14 0.11 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------- ------- --------- -------------- --------- --------------
No psychotic disorder 13014 99.72 99.72 99.37 99.37
Psychotic disorder 37 0.28 100.00 0.28 99.66
<NA> 45 0.34 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------- ------ --------- -------------- --------- --------------
No psychotic disorder 6922 97.85 97.85 80.01 80.01
Psychotic disorder 152 2.15 100.00 1.76 81.77
<NA> 1577 18.23 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$mhd.mania_hypomania_bipolar_or_manicdepression
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------------------------------- ------ --------- -------------- --------- --------------
Not Mania, hypomania, bipolar or manic-depression 62 95.38 95.38 95.38 95.38
Mania, hypomania, bipolar or manic-depression 3 4.62 100.00 4.62 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Mania, hypomania, bipolar or manic-depression 11778 93.09 93.09 93.08 93.08
Mania, hypomania, bipolar or manic-depression 874 6.91 100.00 6.91 99.99
<NA> 1 0.01 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Mania, hypomania, bipolar or manic-depression 12996 99.45 99.45 99.24 99.24
Mania, hypomania, bipolar or manic-depression 72 0.55 100.00 0.55 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------------------------------- ------ --------- -------------- --------- --------------
Not Mania, hypomania, bipolar or manic-depression 7320 96.83 96.83 84.61 84.61
Mania, hypomania, bipolar or manic-depression 240 3.17 100.00 2.77 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$mhd.posttraumatic_stress_disorder_ptsd
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------------------------- ------ --------- -------------- --------- --------------
Not Post-traumatic stress disorder (PTSD) 59 90.77 90.77 90.77 90.77
Post-traumatic stress disorder (PTSD) 6 9.23 100.00 9.23 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------------------------- ------- --------- -------------- --------- --------------
Not Post-traumatic stress disorder (PTSD) 10801 85.37 85.37 85.36 85.36
Post-traumatic stress disorder (PTSD) 1851 14.63 100.00 14.63 99.99
<NA> 1 0.01 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------------------------- ------- --------- -------------- --------- --------------
Not Post-traumatic stress disorder (PTSD) 12837 98.23 98.23 98.02 98.02
Post-traumatic stress disorder (PTSD) 231 1.77 100.00 1.76 99.79
<NA> 28 0.21 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------------------------------- ------ --------- -------------- --------- --------------
Not Post-traumatic stress disorder (PTSD) 6859 90.73 90.73 79.29 79.29
Post-traumatic stress disorder (PTSD) 701 9.27 100.00 8.10 87.39
<NA> 1091 12.61 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$autism_spectrum_disorder
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------- ------ --------- -------------- --------- --------------
No autism spectrum disorder 59 90.77 90.77 90.77 90.77
Autism spectrum disorder 6 9.23 100.00 9.23 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------- ------- --------- -------------- --------- --------------
No autism spectrum disorder 12255 96.96 96.96 96.85 96.85
Autism spectrum disorder 384 3.04 100.00 3.03 99.89
<NA> 14 0.11 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------- ------- --------- -------------- --------- --------------
No autism spectrum disorder 13005 99.65 99.65 99.31 99.31
Autism spectrum disorder 46 0.35 100.00 0.35 99.66
<NA> 45 0.34 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
--------------------------------- ------ --------- -------------- --------- --------------
No autism spectrum disorder 6883 97.30 97.30 79.56 79.56
Autism spectrum disorder 191 2.70 100.00 2.21 81.77
<NA> 1577 18.23 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$mhd.attention_deficit_hyperactivity_disorder
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------------------------------------------------------- ------ --------- -------------- --------- --------------
Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 63 96.92 96.92 96.92 96.92
Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 2 3.08 100.00 3.08 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 12438 98.41 98.41 98.30 98.30
Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 201 1.59 100.00 1.59 99.89
<NA> 14 0.11 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------------------------------------------------------- ------- --------- -------------- --------- --------------
Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 13019 99.75 99.75 99.41 99.41
Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 32 0.25 100.00 0.24 99.66
<NA> 45 0.34 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
-------------------------------------------------------------------------------------- ------ --------- -------------- --------- --------------
Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 6960 98.39 98.39 80.45 80.45
Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD) 114 1.61 100.00 1.32 81.77
<NA> 1577 18.23 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$mhd.personality_disorder
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------ ------ --------- -------------- --------- --------------
Not Personality disorder 59 90.77 90.77 90.77 90.77
Personality disorder 6 9.23 100.00 9.23 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------ ------- --------- -------------- --------- --------------
Not Personality disorder 11688 92.48 92.48 92.37 92.37
Personality disorder 951 7.52 100.00 7.52 99.89
<NA> 14 0.11 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------ ------- --------- -------------- --------- --------------
Not Personality disorder 12975 99.42 99.42 99.08 99.08
Personality disorder 76 0.58 100.00 0.58 99.66
<NA> 45 0.34 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
------------------------------ ------ --------- -------------- --------- --------------
Not Personality disorder 6779 95.83 95.83 78.36 78.36
Personality disorder 295 4.17 100.00 3.41 81.77
<NA> 1577 18.23 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$control
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------ --------- -------------- --------- --------------
No control 64 98.46 98.46 98.46 98.46
Control 1 1.54 100.00 1.54 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
No control 12485 98.90 98.90 98.67 98.67
Control 139 1.10 100.00 1.10 99.77
<NA> 29 0.23 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------- --------- -------------- --------- --------------
No control 3812 29.19 29.19 29.11 29.11
Control 9247 70.81 100.00 70.61 99.72
<NA> 37 0.28 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
---------------- ------ --------- -------------- --------- --------------
No control 4987 66.27 66.27 57.65 57.65
Control 2538 33.73 100.00 29.34 86.98
<NA> 1126 13.02 100.00
Total 8651 100.00 100.00 100.00 100.00
Frequencies
dat$DidNotRespondMHD
Type: Factor
Group: Sample = EDGI
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------- ------ --------- -------------- --------- --------------
Prefer not to say 0 0.00 0.00 0.00 0.00
Responded to MHD 65 100.00 100.00 100.00 100.00
<NA> 0 0.00 100.00
Total 65 100.00 100.00 100.00 100.00
Group: Sample = GLAD
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------- ------- --------- -------------- --------- --------------
Prefer not to say 1 0.01 0.01 0.01 0.01
Responded to MHD 12652 99.99 100.00 99.99 100.00
<NA> 0 0.00 100.00
Total 12653 100.00 100.00 100.00 100.00
Group: Sample = NBR
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------- ------- --------- -------------- --------- --------------
Prefer not to say 0 0.00 0.00 0.00 0.00
Responded to MHD 13096 100.00 100.00 100.00 100.00
<NA> 0 0.00 100.00
Total 13096 100.00 100.00 100.00 100.00
Group: Sample = RAMP
Freq % Valid % Valid Cum. % Total % Total Cum.
----------------------- ------ --------- -------------- --------- --------------
Prefer not to say 1085 12.54 12.54 12.54 12.54
Responded to MHD 7566 87.46 100.00 87.46 100.00
<NA> 0 0.00 100.00
Total 8651 100.00 100.00 100.00 100.00
++MD: These ‘data availability’ chunks caused my R to crash, so I’m hashing them out for now. Others can test if it works for them and unhash as needed, otherwise I think they should be removed.
Data availability GAD
Data availability PHQ
Data availability PCL
Data availability OCI-R
Eating Disorders Genetics Initative Genetic Links to Anxiety and Depression National BioResource Rapid Assessment of MP
Frequencies
dat$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.19 0.19 0.19 0.19
GLAD 12653 36.71 36.90 36.71 36.90
NBR 13096 38.00 74.90 38.00 74.90
RAMP 8651 25.10 100.00 25.10 100.00
<NA> 0 0.00 100.00
Total 34465 100.00 100.00 100.00 100.00
Frequencies
dat$disorders_total_count
Type: Numeric
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
0 12386 38.85 38.85 35.94 35.94
1 5527 17.33 56.18 16.04 51.97
2 8340 26.16 82.34 24.20 76.17
3 3337 10.47 92.81 9.68 85.86
4 1393 4.37 97.17 4.04 89.90
5 569 1.78 98.96 1.65 91.55
6 239 0.75 99.71 0.69 92.24
7 64 0.20 99.91 0.19 92.43
8 23 0.07 99.98 0.07 92.49
9 6 0.02 100.00 0.02 92.51
<NA> 2581 7.49 100.00
Total 34465 100.00 100.00 100.00 100.00
#HR
#SB - Impute no change prepandemic scores
[1] 21 11 16 10 1 4 13 0 14 11 0 9 16 8 7 14 14 6 4 5 9 9 21 19 4
[26] 6 7 9 7 14
[1] NA 2 9 5 NA 21 NA NA NA 3 NA 3 NA 3 3 13 NA NA 3 0 NA 5 15 17 NA
[26] NA 2 NA 4 4
[1] No different Much worse A little worse A little worse No different
[6] Much better No different No different No different A little worse
5 Levels: Much worse A little worse No different ... Much better
[1] 3 1 2 2 3 5 3 3 3 2
[1] 21 11 16 10 1 4 13 0 14 11 0 9 16 8 7 14 14 6 4 5 9 9 21 19 4
[26] 6 7 9 7 14
[1] 21 2 9 5 1 21 13 0 14 3 0 3 16 3 3 13 14 6 3 0 9 5 15 17 4
[26] 6 2 9 4 4
[1] 22 9 15 11 0 9 23 5 23 11 1 20 18 15 4 11 15 19 12 10 17 11 24 27 6
[26] 9 3 11 3 8
[1] 23 1 7 7 NA 20 18 NA 23 5 NA 3 NA 7 3 11 NA NA 12 2 NA 6 19 27 NA
[26] NA 1 NA NA 3
[1] A little worse Much worse A little worse A little worse
[5] No different A little better A little worse No different
[9] A little worse A little worse
5 Levels: Much worse A little worse No different ... Much better
[1] 2 1 2 2 3 4 2 3 2 2
[1] 22 9 15 11 0 9 23 5 23 11 1 20 18 15 4 11 15 19 12 10 17 11 24 27 6
[26] 9 3 11 3 8
[1] 23 1 7 7 0 20 18 5 23 5 1 3 18 7 3 11 15 19 12 2 17 6 19 27 6
[26] 9 1 11 3 3
#Baseline survey completion (N per month)
[1] "2020-04-30 20:51:40 UTC" "2020-04-30 18:55:12 UTC"
[3] "2020-04-30 19:02:04 UTC" "2020-04-30 18:59:37 UTC"
[5] "2020-05-01 10:02:34 UTC" "2020-04-30 22:16:01 UTC"
[7] "2020-04-30 20:36:37 UTC" "2020-04-30 19:02:03 UTC"
[9] "2020-04-30 18:57:33 UTC" "2020-05-01 08:54:01 UTC"
[11] "2020-04-30 21:47:03 UTC" "2020-04-30 19:10:42 UTC"
[13] "2020-04-30 19:29:36 UTC" "2020-04-30 19:20:42 UTC"
[15] "2020-05-01 15:09:23 UTC" "2020-04-30 18:59:55 UTC"
[17] "2020-05-04 16:40:21 UTC" "2020-05-01 16:06:45 UTC"
[19] "2020-05-19 10:34:16 UTC" "2020-04-30 19:11:26 UTC"
Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999
Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999
Warning in hetcor.data.frame(as.data.frame(dat[, mhd.items]), use =
"pairwise.complete.obs"): the correlation matrix has been adjusted to make it
positive-definite
Two-Step Estimates
Correlations/Type of Correlation:
depressive_disorders
depressive_disorders 1
anxiety_disorders 0.8377
eating_disorders 0.5343
obsessive_compulsive_disorders 0.5269
psychotic_disorders 0.5033
mhd.mania_hypomania_bipolar_or_manicdepression 0.5367
mhd.posttraumatic_stress_disorder_ptsd 0.5807
autism_spectrum_disorder 0.3882
mhd.attention_deficit_hyperactivity_disorder 0.3696
mhd.personality_disorder 0.6383
control -0.8096
anxiety_disorders
depressive_disorders Polychoric
anxiety_disorders 1
eating_disorders 0.4969
obsessive_compulsive_disorders 0.6313
psychotic_disorders 0.4568
mhd.mania_hypomania_bipolar_or_manicdepression 0.431
mhd.posttraumatic_stress_disorder_ptsd 0.5833
autism_spectrum_disorder 0.4385
mhd.attention_deficit_hyperactivity_disorder 0.4147
mhd.personality_disorder 0.5838
control -0.8075
eating_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders 1
obsessive_compulsive_disorders 0.5707
psychotic_disorders 0.3641
mhd.mania_hypomania_bipolar_or_manicdepression 0.3653
mhd.posttraumatic_stress_disorder_ptsd 0.3867
autism_spectrum_disorder 0.2962
mhd.attention_deficit_hyperactivity_disorder 0.3309
mhd.personality_disorder 0.5191
control -0.6858
obsessive_compulsive_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders 1
psychotic_disorders 0.3917
mhd.mania_hypomania_bipolar_or_manicdepression 0.3917
mhd.posttraumatic_stress_disorder_ptsd 0.4481
autism_spectrum_disorder 0.3976
mhd.attention_deficit_hyperactivity_disorder 0.3908
mhd.personality_disorder 0.502
control -0.7346
psychotic_disorders
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders 1
mhd.mania_hypomania_bipolar_or_manicdepression 0.7086
mhd.posttraumatic_stress_disorder_ptsd 0.4476
autism_spectrum_disorder 0.2901
mhd.attention_deficit_hyperactivity_disorder 0.3375
mhd.personality_disorder 0.6044
control -0.701
mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression 1
mhd.posttraumatic_stress_disorder_ptsd 0.3492
autism_spectrum_disorder 0.2813
mhd.attention_deficit_hyperactivity_disorder 0.342
mhd.personality_disorder 0.5855
control -0.6882
mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd 1
autism_spectrum_disorder 0.2916
mhd.attention_deficit_hyperactivity_disorder 0.2851
mhd.personality_disorder 0.5077
control -0.6867
autism_spectrum_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder 1
mhd.attention_deficit_hyperactivity_disorder 0.5726
mhd.personality_disorder 0.3777
control -0.6313
mhd.attention_deficit_hyperactivity_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder 1
mhd.personality_disorder 0.3698
control -0.6354
mhd.personality_disorder
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder Polychoric
mhd.personality_disorder 1
control -0.7725
control
depressive_disorders Polychoric
anxiety_disorders Polychoric
eating_disorders Polychoric
obsessive_compulsive_disorders Polychoric
psychotic_disorders Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd Polychoric
autism_spectrum_disorder Polychoric
mhd.attention_deficit_hyperactivity_disorder Polychoric
mhd.personality_disorder Polychoric
control 1
Standard Errors/Numbers of Observations:
depressive_disorders
depressive_disorders 33345
anxiety_disorders 0.004037
eating_disorders 0.01226
obsessive_compulsive_disorders 0.0117
psychotic_disorders 0.02142
mhd.mania_hypomania_bipolar_or_manicdepression 0.01569
mhd.posttraumatic_stress_disorder_ptsd 0.01078
autism_spectrum_disorder 0.02049
mhd.attention_deficit_hyperactivity_disorder 0.02606
mhd.personality_disorder 0.01498
control Inf
anxiety_disorders
depressive_disorders 32428
anxiety_disorders 32428
eating_disorders 0.01209
obsessive_compulsive_disorders 0.01054
psychotic_disorders 0.0202
mhd.mania_hypomania_bipolar_or_manicdepression 0.01546
mhd.posttraumatic_stress_disorder_ptsd 0.01032
autism_spectrum_disorder 0.01985
mhd.attention_deficit_hyperactivity_disorder 0.02525
mhd.personality_disorder 0.01397
control Inf
eating_disorders
depressive_disorders 32836
anxiety_disorders 31919
eating_disorders 32870
obsessive_compulsive_disorders 0.01299
psychotic_disorders 0.02372
mhd.mania_hypomania_bipolar_or_manicdepression 0.01942
mhd.posttraumatic_stress_disorder_ptsd 0.01538
autism_spectrum_disorder 0.02513
mhd.attention_deficit_hyperactivity_disorder 0.02924
mhd.personality_disorder 0.0157
control 1.851
obsessive_compulsive_disorders
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 0.02248
mhd.mania_hypomania_bipolar_or_manicdepression 0.01825
mhd.posttraumatic_stress_disorder_ptsd 0.01397
autism_spectrum_disorder 0.02216
mhd.attention_deficit_hyperactivity_disorder 0.0268
mhd.personality_disorder 0.01562
control 11.73
psychotic_disorders
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 0.01621
mhd.posttraumatic_stress_disorder_ptsd 0.02076
autism_spectrum_disorder 0.03517
mhd.attention_deficit_hyperactivity_disorder 0.03917
mhd.personality_disorder 0.01887
control 3.285
mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 32804
mhd.mania_hypomania_bipolar_or_manicdepression 33345
mhd.posttraumatic_stress_disorder_ptsd 0.01847
autism_spectrum_disorder 0.02949
mhd.attention_deficit_hyperactivity_disorder 0.03287
mhd.personality_disorder 0.01632
control 1.964
mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders 33345
anxiety_disorders 32428
eating_disorders 32836
obsessive_compulsive_disorders 33345
psychotic_disorders 32804
mhd.mania_hypomania_bipolar_or_manicdepression 33345
mhd.posttraumatic_stress_disorder_ptsd 33345
autism_spectrum_disorder 0.02388
mhd.attention_deficit_hyperactivity_disorder 0.02903
mhd.personality_disorder 0.01514
control 1.983
autism_spectrum_disorder
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 32804
mhd.posttraumatic_stress_disorder_ptsd 32804
autism_spectrum_disorder 32829
mhd.attention_deficit_hyperactivity_disorder 0.02798
mhd.personality_disorder 0.02533
control 3.353
mhd.attention_deficit_hyperactivity_disorder
depressive_disorders 32804
anxiety_disorders 31887
eating_disorders 32826
obsessive_compulsive_disorders 32804
psychotic_disorders 32829
mhd.mania_hypomania_bipolar_or_manicdepression 32804
mhd.posttraumatic_stress_disorder_ptsd 32804
autism_spectrum_disorder 32829
mhd.attention_deficit_hyperactivity_disorder 32829
mhd.personality_disorder 0.03046
control 2.038
mhd.personality_disorder control
depressive_disorders 32804 33271
anxiety_disorders 31887 32382
eating_disorders 32826 32800
obsessive_compulsive_disorders 32804 33271
psychotic_disorders 32829 32777
mhd.mania_hypomania_bipolar_or_manicdepression 32804 33271
mhd.posttraumatic_stress_disorder_ptsd 32804 33271
autism_spectrum_disorder 32829 32777
mhd.attention_deficit_hyperactivity_disorder 32829 32777
mhd.personality_disorder 32829 32777
control 1.749 33273
Create the correlation matrix plot using the corrplot package
#Correlations of retro impute and prepan ##Create function to use for all of the corplots
###Whole sample prepan retro correlation
Frequencies
prepan_retro_dat$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.5143 1
phq.sum_score_prepan 0.6904 0.4485
phq.sum_score_retro 0.4443 0.7445
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5886 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 12692 9063
gad.sum_score_retro 0.007727 9082
phq.sum_score_prepan 0.004652 0.008394
phq.sum_score_retro 0.008274 0.004961
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 12662 9411
gad.sum_score_retro 9058 8074
phq.sum_score_prepan 12681 9405
phq.sum_score_retro 0.006739 9431
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 6.754e-309
phq.sum_score_prepan 2.342e-223 2.296e-275
phq.sum_score_retro 4.022e-131 7.718e-221
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 1.101e-125
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
whole_retro_prepan_scatter$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 33 rows containing non-finite values (stat_smooth).
Warning: Removed 33 rows containing missing values (geom_point).
Warning: Removed 33 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 42 rows containing non-finite values (stat_smooth).
Warning: Removed 42 rows containing missing values (geom_point).
Warning: Removed 42 rows containing missing values (geom_point).
###Little/much worse prepan retro correlation
[1] "gad.sum_score_prepan" "gad.sum_score_retro" "phq.sum_score_prepan"
[4] "phq.sum_score_retro"
Frequencies
prepan_retro_dat_worse$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
EDGI 42 0.54 0.54 0.54 0.54
GLAD 7795 99.46 100.00 99.46 100.00
<NA> 0 0.00 100.00
Total 7837 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.5476 1
phq.sum_score_prepan 0.69 0.4785
phq.sum_score_retro 0.4592 0.7366
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5976 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 7820 7808
gad.sum_score_retro 0.007924 7825
phq.sum_score_prepan 0.005931 0.008729
phq.sum_score_retro 0.009446 0.005473
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 7805 6982
gad.sum_score_retro 7805 6988
phq.sum_score_prepan 7817 6979
phq.sum_score_retro 0.007697 6997
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 8.019e-142
phq.sum_score_prepan 1.022e-144 2.245e-111
phq.sum_score_retro 1.966e-123 1.629e-103
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 2.61e-134
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
little_much_worse_retro_prepan$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 14 rows containing non-finite values (stat_smooth).
Warning: Removed 14 rows containing missing values (geom_point).
Warning: Removed 14 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 15 rows containing non-finite values (stat_smooth).
Warning: Removed 15 rows containing missing values (geom_point).
Warning: Removed 15 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 29 rows containing non-finite values (stat_smooth).
Warning: Removed 29 rows containing missing values (geom_point).
Warning: Removed 29 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 11 rows containing non-finite values (stat_smooth).
Warning: Removed 11 rows containing missing values (geom_point).
Warning: Removed 11 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 21 rows containing non-finite values (stat_smooth).
Warning: Removed 21 rows containing missing values (geom_point).
Warning: Removed 21 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 32 rows containing non-finite values (stat_smooth).
Warning: Removed 32 rows containing missing values (geom_point).
Warning: Removed 32 rows containing missing values (geom_point).
###Little/much better prepan retro correlation
[1] "gad.sum_score_prepan" "gad.sum_score_retro" "phq.sum_score_prepan"
[4] "phq.sum_score_retro"
Frequencies
prepan_retro_dat_better$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------ --------- -------------- --------- --------------
EDGI 7 0.56 0.56 0.56 0.56
GLAD 1252 99.44 100.00 99.44 100.00
<NA> 0 0.00 100.00
Total 1259 100.00 100.00 100.00 100.00
Two-Step Estimates
Correlations/Type of Correlation:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1 Pearson
gad.sum_score_retro 0.4858 1
phq.sum_score_prepan 0.6512 0.4084
phq.sum_score_retro 0.3567 0.6977
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan Pearson Pearson
gad.sum_score_retro Pearson Pearson
phq.sum_score_prepan 1 Pearson
phq.sum_score_retro 0.5346 1
Standard Errors/Numbers of Observations:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan 1257 1255
gad.sum_score_retro 0.02158 1257
phq.sum_score_prepan 0.01628 0.02355
phq.sum_score_retro 0.02651 0.01559
phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan 1253 1085
gad.sum_score_retro 1253 1086
phq.sum_score_prepan 1255 1084
phq.sum_score_retro 0.02171 1087
P-values for Tests of Bivariate Normality:
gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan
gad.sum_score_retro 2.716e-29
phq.sum_score_prepan 7.975e-34 3.25e-08
phq.sum_score_retro 1.024e-27 1.832e-07
phq.sum_score_prepan
gad.sum_score_prepan
gad.sum_score_retro
phq.sum_score_prepan
phq.sum_score_retro 1.275e-12
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, :
Finding confidence intervals from a correlation matrix, but n is not specified, arbitrarily set to 1000
Frequencies
little_much_better_retro_prepan$Sample
Type: Character
Freq % Valid % Valid Cum. % Total % Total Cum.
----------- ------- --------- -------------- --------- --------------
EDGI 65 0.51 0.51 0.51 0.51
GLAD 12653 99.49 100.00 99.49 100.00
<NA> 0 0.00 100.00
Total 12718 100.00 100.00 100.00 100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 3 rows containing non-finite values (stat_smooth).
Warning: Removed 3 rows containing missing values (geom_point).
Warning: Removed 3 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 1 rows containing non-finite values (stat_smooth).
Warning: Removed 1 rows containing missing values (geom_point).
Warning: Removed 1 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 4 rows containing non-finite values (stat_smooth).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 4 rows containing non-finite values (stat_smooth).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Removed 4 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 6 rows containing non-finite values (stat_smooth).
Warning: Removed 6 rows containing missing values (geom_point).
Warning: Removed 6 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 10 rows containing non-finite values (stat_smooth).
Warning: Removed 10 rows containing missing values (geom_point).
Warning: Removed 10 rows containing missing values (geom_point).
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
100% bar per diagnosis with percentages of contribution per sample
| Disorder_inc | Sample | n | Prop |
|---|---|---|---|
| Anxiety disorders | EDGI | 38 | 0.00 |
| Anxiety disorders | GLAD | 10098 | 0.66 |
| Anxiety disorders | NBR | 1765 | 0.12 |
| Anxiety disorders | RAMP | 3350 | 0.22 |
| Depressive disorders | EDGI | 47 | 0.00 |
| Depressive disorders | GLAD | 11455 | 0.63 |
| Depressive disorders | NBR | 2756 | 0.15 |
| Depressive disorders | RAMP | 3826 | 0.21 |
| Eating disorders | EDGI | 61 | 0.03 |
| Eating disorders | GLAD | 1384 | 0.65 |
| Eating disorders | NBR | 195 | 0.09 |
| Eating disorders | RAMP | 482 | 0.23 |
| OCDs | EDGI | 18 | 0.01 |
| OCDs | GLAD | 1627 | 0.68 |
| OCDs | NBR | 145 | 0.06 |
| OCDs | RAMP | 607 | 0.25 |
| Psychotic disorders | GLAD | 401 | 0.68 |
| Psychotic disorders | NBR | 37 | 0.06 |
| Psychotic disorders | RAMP | 152 | 0.26 |
| Bipolar disorders | EDGI | 3 | 0.00 |
| Bipolar disorders | GLAD | 874 | 0.74 |
| Bipolar disorders | NBR | 72 | 0.06 |
| Bipolar disorders | RAMP | 240 | 0.20 |
| PTSD | EDGI | 6 | 0.00 |
| PTSD | GLAD | 1851 | 0.66 |
| PTSD | NBR | 231 | 0.08 |
| PTSD | RAMP | 701 | 0.25 |
| Autism | EDGI | 6 | 0.01 |
| Autism | GLAD | 384 | 0.61 |
| Autism | NBR | 46 | 0.07 |
| Autism | RAMP | 191 | 0.30 |
| ADHD | EDGI | 2 | 0.01 |
| ADHD | GLAD | 201 | 0.58 |
| ADHD | NBR | 32 | 0.09 |
| ADHD | RAMP | 114 | 0.33 |
| Personality disorders | EDGI | 6 | 0.00 |
| Personality disorders | GLAD | 951 | 0.72 |
| Personality disorders | NBR | 76 | 0.06 |
| Personality disorders | RAMP | 295 | 0.22 |
| Controls | EDGI | 1 | 0.00 |
| Controls | GLAD | 139 | 0.01 |
| Controls | NBR | 9247 | 0.78 |
| Controls | RAMP | 2538 | 0.21 |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
Age bar plot for whole sample (not grouped) Summary dataset for age and display it as a bar plot
| age_category | n | Prop | order |
|---|---|---|---|
| 16-18 | 1055 | 0.03 | 1 |
| 19-25 | 2617 | 0.08 | 2 |
| 26-35 | 5112 | 0.15 | 3 |
| 36-45 | 4781 | 0.14 | 4 |
| 46-55 | 6988 | 0.20 | 5 |
| 56-65 | 8015 | 0.23 | 6 |
| 66-70 | 3189 | 0.09 | 7 |
| 71-75 | 2010 | 0.06 | 8 |
| 76-80 | 544 | 0.02 | 9 |
| 81-85 | 116 | 0.00 | 10 |
| 86-90 | 27 | 0.00 | 11 |
| 91-100 | 6 | 0.00 | 12 |
| 100+ | 5 | 0.00 | 13 |
summary dataset with frequencies
| highest_education | n | freq |
|---|---|---|
| GCSE/CSE | 4224 | 0.12 |
| NVQ | 2148 | 0.06 |
| A-levels | 6742 | 0.20 |
| University | 19948 | 0.58 |
pie chart labels
Warning: Ignoring unknown parameters: face
Warning: Ignoring unknown parameters: face
Warning: Ignoring unknown aesthetics: face
pie chart - no labels
| Ethnicity_collapsed | n | freq |
|---|---|---|
| European | 30201 | 0.96 |
| Mixed or multiple ethnic origins | 571 | 0.02 |
| Asian or Asian British | 394 | 0.01 |
| African or African British | 159 | 0.01 |
| Other | 193 | 0.01 |
pie chart with only European label
Warning: Ignoring unknown parameters: face
Pie chart - no labels
| Gender_collapsed | n | freq |
|---|---|---|
| Male | 9698 | 0.28 |
| Female | 24387 | 0.71 |
| Non-binary/Self-defined | 380 | 0.01 |
Warning: Ignoring unknown parameters: face
Pie chart - only labelling two largest groups (male and female)
Warning: Ignoring unknown parameters: face
Pie chart - no labels
+++CH: Diagnosis, ethnicity, gender, age
| age_category | Sample | n | Prop |
|---|---|---|---|
| 16-18 | EDGI | 7 | 0.00 |
| 16-18 | GLAD | 341 | 0.01 |
| 16-18 | NBR | 10 | 0.00 |
| 16-18 | RAMP | 697 | 0.02 |
| 19-25 | EDGI | 17 | 0.00 |
| 19-25 | GLAD | 1537 | 0.04 |
| 19-25 | NBR | 291 | 0.01 |
| 19-25 | RAMP | 772 | 0.02 |
| 26-35 | EDGI | 21 | 0.00 |
| 26-35 | GLAD | 2850 | 0.08 |
| 26-35 | NBR | 1165 | 0.03 |
| 26-35 | RAMP | 1076 | 0.03 |
| 36-45 | EDGI | 8 | 0.00 |
| 36-45 | GLAD | 2459 | 0.07 |
| 36-45 | NBR | 1440 | 0.04 |
| 36-45 | RAMP | 874 | 0.03 |
| 46-55 | EDGI | 4 | 0.00 |
| 46-55 | GLAD | 2821 | 0.08 |
| 46-55 | NBR | 2644 | 0.08 |
| 46-55 | RAMP | 1519 | 0.04 |
| 56-65 | EDGI | 5 | 0.00 |
| 56-65 | GLAD | 1908 | 0.06 |
| 56-65 | NBR | 3901 | 0.11 |
| 56-65 | RAMP | 2201 | 0.06 |
| 66-70 | EDGI | 3 | 0.00 |
| 66-70 | GLAD | 433 | 0.01 |
| 66-70 | NBR | 1967 | 0.06 |
| 66-70 | RAMP | 786 | 0.02 |
| 71-75 | GLAD | 231 | 0.01 |
| 71-75 | NBR | 1265 | 0.04 |
| 71-75 | RAMP | 514 | 0.01 |
| 76-80 | GLAD | 59 | 0.00 |
| 76-80 | NBR | 334 | 0.01 |
| 76-80 | RAMP | 151 | 0.00 |
| 81-85 | GLAD | 13 | 0.00 |
| 81-85 | NBR | 61 | 0.00 |
| 81-85 | RAMP | 42 | 0.00 |
| 86-90 | GLAD | 1 | 0.00 |
| 86-90 | NBR | 13 | 0.00 |
| 86-90 | RAMP | 13 | 0.00 |
| 91-100 | NBR | 2 | 0.00 |
| 91-100 | RAMP | 4 | 0.00 |
| 100+ | NBR | 3 | 0.00 |
| 100+ | RAMP | 2 | 0.00 |
| Sample | age_category | n | Prop |
|---|---|---|---|
| EDGI | 16-18 | 7 | 0.11 |
| EDGI | 19-25 | 17 | 0.26 |
| EDGI | 26-35 | 21 | 0.32 |
| EDGI | 36-45 | 8 | 0.12 |
| EDGI | 46-55 | 4 | 0.06 |
| EDGI | 56-65 | 5 | 0.08 |
| EDGI | 66-70 | 3 | 0.05 |
| GLAD | 16-18 | 341 | 0.03 |
| GLAD | 19-25 | 1537 | 0.12 |
| GLAD | 26-35 | 2850 | 0.23 |
| GLAD | 36-45 | 2459 | 0.19 |
| GLAD | 46-55 | 2821 | 0.22 |
| GLAD | 56-65 | 1908 | 0.15 |
| GLAD | 66-70 | 433 | 0.03 |
| GLAD | 71-75 | 231 | 0.02 |
| GLAD | 76-80 | 59 | 0.00 |
| GLAD | 81-85 | 13 | 0.00 |
| GLAD | 86-90 | 1 | 0.00 |
| NBR | 16-18 | 10 | 0.00 |
| NBR | 19-25 | 291 | 0.02 |
| NBR | 26-35 | 1165 | 0.09 |
| NBR | 36-45 | 1440 | 0.11 |
| NBR | 46-55 | 2644 | 0.20 |
| NBR | 56-65 | 3901 | 0.30 |
| NBR | 66-70 | 1967 | 0.15 |
| NBR | 71-75 | 1265 | 0.10 |
| NBR | 76-80 | 334 | 0.03 |
| NBR | 81-85 | 61 | 0.00 |
| NBR | 86-90 | 13 | 0.00 |
| NBR | 91-100 | 2 | 0.00 |
| NBR | 100+ | 3 | 0.00 |
| RAMP | 16-18 | 697 | 0.08 |
| RAMP | 19-25 | 772 | 0.09 |
| RAMP | 26-35 | 1076 | 0.12 |
| RAMP | 36-45 | 874 | 0.10 |
| RAMP | 46-55 | 1519 | 0.18 |
| RAMP | 56-65 | 2201 | 0.25 |
| RAMP | 66-70 | 786 | 0.09 |
| RAMP | 71-75 | 514 | 0.06 |
| RAMP | 76-80 | 151 | 0.02 |
| RAMP | 81-85 | 42 | 0.00 |
| RAMP | 86-90 | 13 | 0.00 |
| RAMP | 91-100 | 4 | 0.00 |
| RAMP | 100+ | 2 | 0.00 |
| age_category | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| 16-18 | 7 | 341 | 10 | 697 | 0.11 | 0.03 | 0.00 | 0.08 |
| 19-25 | 17 | 1537 | 291 | 772 | 0.26 | 0.12 | 0.02 | 0.09 |
| 26-35 | 21 | 2850 | 1165 | 1076 | 0.32 | 0.23 | 0.09 | 0.12 |
| 36-45 | 8 | 2459 | 1440 | 874 | 0.12 | 0.19 | 0.11 | 0.10 |
| 46-55 | 4 | 2821 | 2644 | 1519 | 0.06 | 0.22 | 0.20 | 0.18 |
| 56-65 | 5 | 1908 | 3901 | 2201 | 0.08 | 0.15 | 0.30 | 0.25 |
| 66-70 | 3 | 433 | 1967 | 786 | 0.05 | 0.03 | 0.15 | 0.09 |
| 71-75 | NA | 231 | 1265 | 514 | NA | 0.02 | 0.10 | 0.06 |
| 76-80 | NA | 59 | 334 | 151 | NA | 0.00 | 0.03 | 0.02 |
| 81-85 | NA | 13 | 61 | 42 | NA | 0.00 | 0.00 | 0.00 |
| 86-90 | NA | 1 | 13 | 13 | NA | 0.00 | 0.00 | 0.00 |
| 91-100 | NA | NA | 2 | 4 | NA | NA | 0.00 | 0.00 |
| 100+ | NA | NA | 3 | 2 | NA | NA | 0.00 | 0.00 |
| age_category | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| 16-18 | 7 | 341 | 10 | 697 |
| 19-25 | 17 | 1537 | 291 | 772 |
| 26-35 | 21 | 2850 | 1165 | 1076 |
| 36-45 | 8 | 2459 | 1440 | 874 |
| 46-55 | 4 | 2821 | 2644 | 1519 |
| 56-65 | 5 | 1908 | 3901 | 2201 |
| 66-70 | 3 | 433 | 1967 | 786 |
| 71-75 | 0 | 231 | 1265 | 514 |
| 76-80 | 0 | 59 | 334 | 151 |
| 81-85 | 0 | 13 | 61 | 42 |
| 86-90 | 0 | 1 | 13 | 13 |
| 91-100 | 0 | 0 | 2 | 4 |
| 100+ | 0 | 0 | 3 | 2 |
Warning in chisq.test(.): Chi-squared approximation may be incorrect
| statistic | p.value | parameter | method |
|---|---|---|---|
| 5856.536 | 0 | 36 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
| Ethnicity_collapsed | Sample | n | Prop |
|---|---|---|---|
| European | EDGI | 60 | 0.00 |
| European | GLAD | 12090 | 0.35 |
| European | NBR | 9927 | 0.29 |
| European | RAMP | 8124 | 0.24 |
| Mixed or multiple ethnic origins | EDGI | 3 | 0.00 |
| Mixed or multiple ethnic origins | GLAD | 276 | 0.01 |
| Mixed or multiple ethnic origins | NBR | 118 | 0.00 |
| Mixed or multiple ethnic origins | RAMP | 174 | 0.01 |
| Asian or Asian British | EDGI | 1 | 0.00 |
| Asian or Asian British | GLAD | 98 | 0.00 |
| Asian or Asian British | NBR | 126 | 0.00 |
| Asian or Asian British | RAMP | 169 | 0.00 |
| African or African British | GLAD | 37 | 0.00 |
| African or African British | NBR | 46 | 0.00 |
| African or African British | RAMP | 76 | 0.00 |
| Other | EDGI | 1 | 0.00 |
| Other | GLAD | 122 | 0.00 |
| Other | RAMP | 70 | 0.00 |
| NA | GLAD | 30 | 0.00 |
| NA | NBR | 2879 | 0.08 |
| NA | RAMP | 38 | 0.00 |
| Sample | Ethnicity_collapsed | n | Prop |
|---|---|---|---|
| EDGI | European | 60 | 0.92 |
| EDGI | Mixed or multiple ethnic origins | 3 | 0.05 |
| EDGI | Asian or Asian British | 1 | 0.02 |
| EDGI | Other | 1 | 0.02 |
| GLAD | European | 12090 | 0.96 |
| GLAD | Mixed or multiple ethnic origins | 276 | 0.02 |
| GLAD | Asian or Asian British | 98 | 0.01 |
| GLAD | African or African British | 37 | 0.00 |
| GLAD | Other | 122 | 0.01 |
| GLAD | NA | 30 | 0.00 |
| NBR | European | 9927 | 0.76 |
| NBR | Mixed or multiple ethnic origins | 118 | 0.01 |
| NBR | Asian or Asian British | 126 | 0.01 |
| NBR | African or African British | 46 | 0.00 |
| NBR | NA | 2879 | 0.22 |
| RAMP | European | 8124 | 0.94 |
| RAMP | Mixed or multiple ethnic origins | 174 | 0.02 |
| RAMP | Asian or Asian British | 169 | 0.02 |
| RAMP | African or African British | 76 | 0.01 |
| RAMP | Other | 70 | 0.01 |
| RAMP | NA | 38 | 0.00 |
| Ethnicity_collapsed | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| European | 60 | 12090 | 9927 | 8124 | 0.92 | 0.96 | 0.76 | 0.94 |
| Mixed or multiple ethnic origins | 3 | 276 | 118 | 174 | 0.05 | 0.02 | 0.01 | 0.02 |
| Asian or Asian British | 1 | 98 | 126 | 169 | 0.02 | 0.01 | 0.01 | 0.02 |
| Other | 1 | 122 | NA | 70 | 0.02 | 0.01 | NA | 0.01 |
| African or African British | NA | 37 | 46 | 76 | NA | 0.00 | 0.00 | 0.01 |
| NA | NA | 30 | 2879 | 38 | NA | 0.00 | 0.22 | 0.00 |
| Ethnicity_collapsed | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| European | 60 | 12090 | 9927 | 8124 |
| Mixed or multiple ethnic origins | 3 | 276 | 118 | 174 |
| Asian or Asian British | 1 | 98 | 126 | 169 |
| African or African British | 0 | 37 | 46 | 76 |
| Other | 1 | 122 | 0 | 70 |
| NA | 0 | 30 | 2879 | 38 |
Warning in chisq.test(.): Chi-squared approximation may be incorrect
| statistic | p.value | parameter | method |
|---|---|---|---|
| 5106.32 | 0 | 15 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
| Gender_collapsed | Sample | n | Prop |
|---|---|---|---|
| Male | EDGI | 1 | 0.00 |
| Male | GLAD | 2279 | 0.07 |
| Male | NBR | 5641 | 0.16 |
| Male | RAMP | 1777 | 0.05 |
| Female | EDGI | 62 | 0.00 |
| Female | GLAD | 10127 | 0.29 |
| Female | NBR | 7416 | 0.22 |
| Female | RAMP | 6782 | 0.20 |
| Non-binary/Self-defined | EDGI | 2 | 0.00 |
| Non-binary/Self-defined | GLAD | 247 | 0.01 |
| Non-binary/Self-defined | NBR | 39 | 0.00 |
| Non-binary/Self-defined | RAMP | 92 | 0.00 |
| Sample | Gender_collapsed | n | Prop |
|---|---|---|---|
| EDGI | Male | 1 | 0.02 |
| EDGI | Female | 62 | 0.95 |
| EDGI | Non-binary/Self-defined | 2 | 0.03 |
| GLAD | Male | 2279 | 0.18 |
| GLAD | Female | 10127 | 0.80 |
| GLAD | Non-binary/Self-defined | 247 | 0.02 |
| NBR | Male | 5641 | 0.43 |
| NBR | Female | 7416 | 0.57 |
| NBR | Non-binary/Self-defined | 39 | 0.00 |
| RAMP | Male | 1777 | 0.21 |
| RAMP | Female | 6782 | 0.78 |
| RAMP | Non-binary/Self-defined | 92 | 0.01 |
| Gender_collapsed | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| Male | 1 | 2279 | 5641 | 1777 | 0.02 | 0.18 | 0.43 | 0.21 |
| Female | 62 | 10127 | 7416 | 6782 | 0.95 | 0.80 | 0.57 | 0.78 |
| Non-binary/Self-defined | 2 | 247 | 39 | 92 | 0.03 | 0.02 | 0.00 | 0.01 |
| Gender_collapsed | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| Male | 1 | 2279 | 5641 | 1777 |
| Female | 62 | 10127 | 7416 | 6782 |
| Non-binary/Self-defined | 2 | 247 | 39 | 92 |
Warning in chisq.test(.): Chi-squared approximation may be incorrect
| statistic | p.value | parameter | method |
|---|---|---|---|
| 2455.701 | 0 | 6 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
| highest_education | Sample | n | Prop |
|---|---|---|---|
| GCSE/CSE | EDGI | 7 | 0.00 |
| GCSE/CSE | GLAD | 1078 | 0.03 |
| GCSE/CSE | NBR | 1876 | 0.06 |
| GCSE/CSE | RAMP | 1263 | 0.04 |
| NVQ | EDGI | 3 | 0.00 |
| NVQ | GLAD | 918 | 0.03 |
| NVQ | NBR | 1227 | 0.04 |
| A-levels | EDGI | 20 | 0.00 |
| A-levels | GLAD | 2713 | 0.08 |
| A-levels | NBR | 2320 | 0.07 |
| A-levels | RAMP | 1689 | 0.05 |
| University | EDGI | 35 | 0.00 |
| University | GLAD | 7515 | 0.23 |
| University | NBR | 6983 | 0.21 |
| University | RAMP | 5415 | 0.16 |
| Sample | highest_education | n | Prop |
|---|---|---|---|
| EDGI | GCSE/CSE | 7 | 0.11 |
| EDGI | NVQ | 3 | 0.05 |
| EDGI | A-levels | 20 | 0.31 |
| EDGI | University | 35 | 0.54 |
| GLAD | GCSE/CSE | 1078 | 0.09 |
| GLAD | NVQ | 918 | 0.08 |
| GLAD | A-levels | 2713 | 0.22 |
| GLAD | University | 7515 | 0.61 |
| NBR | GCSE/CSE | 1876 | 0.15 |
| NBR | NVQ | 1227 | 0.10 |
| NBR | A-levels | 2320 | 0.19 |
| NBR | University | 6983 | 0.56 |
| RAMP | GCSE/CSE | 1263 | 0.15 |
| RAMP | A-levels | 1689 | 0.20 |
| RAMP | University | 5415 | 0.65 |
| highest_education | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| GCSE/CSE | 7 | 1078 | 1876 | 1263 | 0.11 | 0.09 | 0.15 | 0.15 |
| NVQ | 3 | 918 | 1227 | NA | 0.05 | 0.08 | 0.10 | NA |
| A-levels | 20 | 2713 | 2320 | 1689 | 0.31 | 0.22 | 0.19 | 0.20 |
| University | 35 | 7515 | 6983 | 5415 | 0.54 | 0.61 | 0.56 | 0.65 |
| highest_education | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| GCSE/CSE | 7 | 1078 | 1876 | 1263 |
| NVQ | 3 | 918 | 1227 | 0 |
| A-levels | 20 | 2713 | 2320 | 1689 |
| University | 35 | 7515 | 6983 | 5415 |
Warning in chisq.test(.): Chi-squared approximation may be incorrect
| statistic | p.value | parameter | method |
|---|---|---|---|
| 1207.855 | 0 | 12 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth
| key_worker | Sample | n | Prop |
|---|---|---|---|
| No key worker | EDGI | 42 | 0.00 |
| No key worker | GLAD | 7249 | 0.21 |
| No key worker | NBR | 8660 | 0.25 |
| No key worker | RAMP | 5841 | 0.17 |
| Key worker | EDGI | 23 | 0.00 |
| Key worker | GLAD | 5366 | 0.16 |
| Key worker | NBR | 4414 | 0.13 |
| Key worker | RAMP | 2765 | 0.08 |
| NA | GLAD | 38 | 0.00 |
| NA | NBR | 22 | 0.00 |
| NA | RAMP | 45 | 0.00 |
| Sample | key_worker | n | Prop |
|---|---|---|---|
| EDGI | No key worker | 42 | 0.65 |
| EDGI | Key worker | 23 | 0.35 |
| GLAD | No key worker | 7249 | 0.57 |
| GLAD | Key worker | 5366 | 0.42 |
| GLAD | NA | 38 | 0.00 |
| NBR | No key worker | 8660 | 0.66 |
| NBR | Key worker | 4414 | 0.34 |
| NBR | NA | 22 | 0.00 |
| RAMP | No key worker | 5841 | 0.68 |
| RAMP | Key worker | 2765 | 0.32 |
| RAMP | NA | 45 | 0.01 |
| key_worker | n_EDGI | n_GLAD | n_NBR | n_RAMP | Prop_EDGI | Prop_GLAD | Prop_NBR | Prop_RAMP |
|---|---|---|---|---|---|---|---|---|
| No key worker | 42 | 7249 | 8660 | 5841 | 0.65 | 0.57 | 0.66 | 0.68 |
| Key worker | 23 | 5366 | 4414 | 2765 | 0.35 | 0.42 | 0.34 | 0.32 |
| NA | NA | 38 | 22 | 45 | NA | 0.00 | 0.00 | 0.01 |
| key_worker | EDGI | GLAD | NBR | RAMP |
|---|---|---|---|---|
| No key worker | 42 | 7249 | 8660 | 5841 |
| Key worker | 23 | 5366 | 4414 | 2765 |
| NA | 0 | 38 | 22 | 45 |
Warning in chisq.test(.): Chi-squared approximation may be incorrect
| statistic | p.value | parameter | method |
|---|---|---|---|
| 332.6547 | 0 | 6 | Pearson’s Chi-squared test |
Warning: Ignoring unknown parameters: binwidth
Warning: Ignoring unknown parameters: binwidth